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Abstract: We consider the superpotential formalism to describe the evolution of D scalar 
fields during inflation, generalizing it to include the case with non-canonical kinetic terms. 
We provide a characterization of the attractor behaviour of the background evolution in 
terms of first and second slow-roll parameters (which need not be small). We find that the 
superpotential is useful in justifying the separate universe approximation from the gradient 
expansion, and also in computing the spectra of primordial perturbations around attractor 
solutions in the 5N formalism. As an application, we consider a class of models where the 
background trajectories for the inflaton helds are derived from a product separable super¬ 
potential. In the perspective of the holographic inflation scenario, such models are dual 
to a deformed CFT boundary theory, with D mutually uncorrelated deformation opera¬ 
tors. We compute the bulk power spectra of primordial adiabatic and entropy cosmological 
perturbations, and show that the results agree with the ones obtained by using conformal 
perturbation theory in the dual picture. 
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1. Introduction 


Measurements of the temperature anisotropies and polarization of the cosmic microwave 
background place stringent constraints on a wide range of inflationary models. While cur¬ 
rent data are consistent with single-field inflation, multi-field scenarios arise quite naturally 
in most attempts to embed inflation within a broader theory, and it is therefore important 
to address this more generic situation. 

The amplitude of primordial perturbations is often described in terms of C, the curva¬ 
ture perturbation on hypersurfaces of constant energy density. In single-field inflation this 
quantity is conserved on super-horizon scales, but in multi-held inhation it can evolve after 
Hubble crossing iil. To take this evolution into account, it is convenient to use the 
so-called 5N formalism iil 0,1 , which gives the correlators of C at the end of inhation 
in terms of the correlators of the scalar-helds huctuations at Hubble crossing. The relation 
between them is established by considering the background evolution of an ensemble of 
homogeneous universes with different initial conditions (for a review, see Refs. i O). 
This approach is not restricted to the correlators of but can also be applied to entropy 
perturbations (see e.g. 0 , 0 ). 

A useful tool for describing the evolution on super-horizon scales is the Hamilton- 
Jacobi (H-J) formalism hrst introduced by Salopek and Bond Q (see also |^). This 
was originally developed for the case of D scalar helds (j)^ (here and below capital latin 
indices run from 1 to D) with canonical kinetic terms, but for generality here 

we shall consider an extension to Lagrangians of the form P{X^^, where 

. ( 1 . 1 ) 


The hrst step is to encode the background dynamics in a time independent H-J equation 
for a “superpotential” W. This object is just the Hubble rate expressed as a function of 
the held values (j)^, 




( 1 . 2 ) 


where a complete solution of the H-J equation contains an equal number of integration 
constants c/, which account for the freedom in the initial values of the held momenta. As 
we shall see, cosmological evolution can then be seen as a “gradient how” of W (0^) in 
held space. The purpose of this paper is to further develop the 6N formalism by taking 
advantage of this description. 

The superpotential approach is also interesting in connection with the possibility of a 
holographic description of inhation [|T^ |l^, |l^. For the case of de Sitter, this idea was 

hrst considered in Refs. |^, ^], by analogy with the gauge/gravity duality which holds in 
asymptotically AdS spaces (see also Refs. |21, 22|). Recently, held theories which are dual 
to de Sitter have been identihed for the case of higher spin gravity On the other hand, 
for the case of Einstein gravity, the duality remains at an exploratory stage. In the absence 
of a more concrete realization, a fruitful strategy has been to focus on small deformations 
of a generic boundary conformal held theory (GET), with Lagrangian of the form L = 
C^CFT + This setup is characterized by a few parameters, such as the central 


- 2 - 





charge of the CFT, and the operator product expansion coefficients for the deformation 
operators O/. Such parametrization allows for some explicit calculations, which can be 

^ . With the identification = (j)^ 


done by using conformal perturbation theory |24, 

(up to a proportionality constant), the renormalisation group (RG) flow of the couplings 
in the boundary theory corresponds to field evolution in the bulk, while the superpotential 
W plays the role of a c-function for the RG flow. 

The case with a single deformation operator O corresponds to single field inflation, 
which has been extensively considered in the literature 

II, |T 


17, \m, 24, 25, 27, 28, 29 


30, 31 


37, 


, m H, H, H|. In particular. Refs. [H, 
studied the power spectrum and bispectrum of C, showing agreement between the boundary 
and bulk calculations.^ The four-point correlation function was also computed in Ref. p3|. 


recovering the result from the bulk calculation of Ref. [^| in the slow-roll regime. 

More recently, Ref. [|6| extended the holographic approach to multi-field inflation by 
considering a GFT with D mutually uncorrelated deformation operators O/, and the pri¬ 
mordial power spectra for adiabatic and entropy perturbations were computed in conformal 
perturbation theory. As an application of the methods presented in this paper, here we 
will compare the results of Ref. with a bulk calculation based on the 5N technique. 

The paper is organized as follows. In Sec. ^ we review the superpotential formalism. 
We also give a characterization of the attractor behavior of background solutions in terms 
of the first and second slow-roll parameters. Note, however, that such parameters will 
not be required to be small in our discussion. In Sec. ^ we review the separate universe 
approximation, on which the 5N formalism is based. As we shall see, the superpotential 
will be very useful for computing the primordial spectra, particularly in the case when the 
background is an attractor. The expressions for the primordial spectra of adiabatic and 
entropy perturbations in terms of the superpotential are given in Sec. |[ In Sec. we 
compare these results with those recently obtained in Ref. from the holographic point 
of view. In Sec. we elaborate on several explicit models of inflation with a product sepa¬ 
rable superpotential, which should be dual to a QFT with uncorrelated multi-deformation 
operators. Finally, we conclude in Sec. H 


2. Superpotential and background evolution 

Consider a scalar field theory in a (d -|- l)-dimensional spacetime, with an action of the 
form 


S = j . (2.1) 

The corresponding energy momentum tensor is given by 

T/.1/ = Pij + Pg^u, (2.2) 

^Furthermore, in Ref. it was shown that the power spectrum of is conserved at large scales, 
as expected in the standard cosmological perturbation theory for the case of single field models. The 
conservation of higher order correlation functions, however, remains an unresolved open issue. 
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where we use the notation 


(2.3) 


PiJ = 


dP 


Since is symmetric about I and J, so is Pjj- If the gradients of all fields are aligned in 
the same time-like direction, the energy momentum tensor (|2.2| ) has the form of a perfect 
fluid, with pressure P and energy density given by 


p = 2PijX^J - P 


(2.4) 


The perfect fluid form will be valid for our background solution, where all fields depend 
only on time. However, in general the fluid will be imperfect for a perturbed solution. 

We assume a flat {d -|- l)-dimensional FRW universe described by the metric ds^ = 
—dt^ -|- a^(t)dx^. The background field equations for (f)'^ {t) are given by 


{Pu^^y+dHPu^^ - {dP/dcl)j) = 0. 


Here 


a 


is the expansion rate. The Friedmann equation reads 


= 


2k? 2k^ 


d{d - 1) ^ d{d - 1) 


{2PijX^J - P) , 


and 




(2.5) 


( 2 . 6 ) 


(2.7) 


( 2 . 8 ) 


1^,11), the 


where = 8ttG is the gravitational constant. 

As pointed out by Salopek and Bond in Ref. |5| (see also Refs, 
field equations for canonical scalar helds can be recast into first-order form by introducing 
a superpotential. Here, we show that the same treatment applies to non-canonical scalar 
fields. Let us start by defining the momenta 


TT/ = 


dP 


= PiJ^' 


(2.9) 


In general Pjj can have explicit dependence in field velocities (f)^, but we assume that our 
Lagrangian is non-singular, so that ( |2.9|) can be solved for the field velocities as a function 
of positions and momenta: 

( 2 . 10 ) 


= fU 


Note that the momenta ttj differ from the canonical ones Hj = dC/dcf)'^ = \/—g nj by 
a factor of but with our choice, the relation o does not contain any explicit time 
dependence through the scale factor a. We may now use (|2.10|) in order to remove all 
occurrences of in the energy density. Then, substituting 


H = 2W{(p 


K\ 


( 2 . 11 ) 
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and 


TTl = Pjjct)'^ = - 


2{d-l)dW 


( 2 . 12 ) 


K? dcf)^ ’ 

which satisfies Eq. ( |2.8D , into the Friedmann equation ( |2.7D , we obtain a partial differential 
equation for the superpotential VE((/)^): 


= 


2d(d-l) 


p [4>^, dW/d<l) 




(2.13) 


Following Ref. |^, we shall call this the separated, or time independent, H-J equation.^ 
The substitutions ( 2.11 ) and (|2.12| ) are actually motivated by the form of the momentum 
constraint in the long wavelength limit ||^. We will discuss this constraint in the next 
Section, where we consider cosmological perturbations [see Eq. ( |3.12| )]. Nonetheless, for the 
time being, we may simply think of ( p.ll ), ( 2.12| ) and (|2.13| ) as a convenient reformulation 
of the background equations of motion. Indeed, in Appendix we show that for any 
solution of (^.13 ), Eqs. ( 2.11 ) and (p.l2 ) generate a solution of the field equations 

(2.5) and O- 

A complete solution of the H-J equation kF(0^,Cif) depending on D independent 
integration constants ck, can be used for generating any background solution from the 
first order equations (2.12) and ( p.llj ). The existence of such a solution can be understood 
as follows. Given some initial data (<^^, vr^) at t = Iq, the background equations of motion 
can be solved in order to find JJ(t; 0^, tt^). If one of the fields has a monotonic evolution, 
we can use it as the time variable, and express the initial positions of the rest in terms of 
their positions at the time t. With this, we have H = H{(j)^] tt^). Finally, we may choose 
ck = invertible relation between ck and tt^), leading to IT = ,ck)/2.. 

2.1 Cosmological evolution as a gradient flow 

A cosmological solution can be thought of as a trajectory in field space, parametrized by the 
e-folding number N. For each field (j)^, we may define the corresponding beta function (3^ 
as the dimensionless component of the tangent vector d/diV in the corresponding direction, 

d(/)^ 


— = , 
dA^ K d(p^ ’ 


(3^ = K 


diV 


(2.14) 


The standard “first” slow-roll parameter can then be written in terms of the beta functions 
as 

H _kHp + P) _ PuP^P^ 


£l = 


> 0 , 


(2.15) 


(d-l)Jj2 {d-1) 

where the last inequality holds provided that the null energy condition is satisfied. In what 
follows, we shall assume that Pjj is non-degenerate. 

In that case, Pij can be thought of as a metric, which we may use in order to raise 
and lower indices. Introducing 

f = -lnW, (2.16) 


^For the case of canonical scalar fields, it is shown in ^ that the separated H-J equation can be obtained 
from the Einstein-Hamilton-Jacobi equation, after factoring out the volume of space in Hamilton’s principal 
function S oc a'^W . 
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Eq. (2.12) can be written as 


(2.17) 


Pi = PiJp^ 


d-ldf 

K ’ 


which expresses the beta functions as a gradient. From Eq. (|2.18) 


dN 


d-1 


£1 > 0 , 


(2.18) 


and therefore / is monotonically increasing along any cosmological trajectory. Here, we 
have used the notation = PiPP Also, since we assume that Pjj is positive definite, the 
last inequality is strict provided that at least one of the fields is moving. 

If we interpret the P^ as beta functions for couplings in the boundary theory, 

then / plays the role of a c-function which decreases monotonically along the RG flow 
(note that the RG flow proceeds from the UV to the IR, which is opposite to the direction 
of cosmological evolution as a function of N). In this interpretation, Pjj corresponds to 
the Zamolodchicov metric in the space of couplings [^. Such metric is needed in order 
to relate the components of the RG flow tangent vector PP which have the upper index, 
to the gradient of a c-function, whose components have the lower index. We shall further 
elaborate on the dual picture in Sec. ^ 


2.2 Attractor behaviour 

As mentioned above, the integration constants ck in W{p^, ck) correspond to the freedom 
of choosing the initial value of the momenta vr^ for given initial values of the fields pQ . 
In the present context, the initial conditions for the long wavelength bulk evolution are 
implemented around the time of horizon crossing. Under a small variation of the integration 
constants, the gradient flow will change. It is then important to characterize whether such 
dependence remains at late times, or whether it decays. 

First of all, it should be noted that along the field trajectory 


dW 

dcK 


A^a-P 


(2.19) 


where are constants and we remind that d is the number of spatial dimensions. For the 
case of scalar fields with canonical kinetic terms, this result was first derived in Ref. [^. We 
show in Appendix ^ that the same result generalizes to Lagrangians with non-canonical 
kinetic terms. Eq. ( 2.1S| ) tells us that the effect of the integration constants on the 
Hubble rate {H = 2W) decays quite rapidly with the scale factor. Of course, this does not 
immediately imply that the gradient flow will always have an attractor behaviour. It is not 
enough that the Hubble rate converges to the unperturbed value, but also the perturbed 
trajectories should converge to the unperturbed ones. It seems difficult to formulate a 
sufficient condition for such convergence in the general case. Here, we shall consider a 
necessary condition, which is also a sufficient condition in the one-field case, or when there 
is effectively just an adiabatic perturbation by the end of inflation. 
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The idea is to consider the change in the held momentum, projected on the unperturbed 
trajectory. The fractional change in under a small variation Ack is given by 

= ( 2 . 20 ) 

Here, we use the notation tt^ = , and (5i and 82 correspond, respectively, to the 

hrst and second terms in the right hand side of the last equality. For the background to 
be an attractor, we require that these go to zero as the universe expands. 


lim (5i 2 = 0. 

a—>-CO ’ 


( 2 . 21 ) 


Since the metric P/j depends on momenta, it will have dependence on ck, and so the two 
terms (5i and <52 can have a somewhat different behaviour. We require that the condition 
(2.21) should hold separately for the two terms, so that the projection of both Avr^ and 
Att/ on the unperturbed trajectory can be considered to be small. 

Using Eq. ( p.l2 ) it is straightforward to show that 

Ack. (2.22) 


/A 4(d-l)^^ d 

, A», =- 




dcK 


Using (2.1£) in (2.22), and (2.7), we have 


tt^Att/ = [dcjiP) Ack = Ap. 


(2.23) 


For later convenience, here we have used the Friedmann equation to write kF^ in terms of 
the energy density p (which depends on the ck only through the kinetic variables 
Using (El) again in order to evaluate the derivative with respect to ck, we have 

(2.24) 


= 


7AW0, 


where AHo = A'^Acj is the change in W at the initial time due to the variation Acj. 

From ( |2.24 ) we find that the first condition (2.21) for the attractor behavior, concerning 
the behavior of (5i, will be satisfied provided that 

d — ei+e2>p>0. (2.25) 

Here we have introduced the second slow-roll parameter 

diV eiP ’ ^ ^ 

and a positive constant p. Eq. (^.25 ) guarantees that (5i decays faster than a~P. If p is 


small, then the approach to the attractor can be slow, requiring a time-scale of the order 
AN ~ p~^ e-foldings. An efficient approach to the attractor on the Hubble time-scale 
requires p > 1. 

In addition, we need to consider the behaviour of 82 , which may place additional 
restrictions. Since = {p + P), from (|2.20|) and (2.23) we immediately obtain 


tt/Att'^ = (dci^P) Ack = AP. 


(2.27) 
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Note that the pressure P only depends on the ck through the kinetic variables X^'\ and 
therefore 

dc^P = Pijdc^X^^. (2.28) 

The same is true of the energy density, 

dc^p = pijdc^X^-^ (2.29) 

where pij = Pij + 2X^^ Pkl,ij is the symmetrized partial derivative of p with respect to 
X^^. Since in general pij ^ Pij, 62 can behave differently than (5i. 

For instance if the background solution satisfies the relation 

AP = cl Ap, (2.30) 


for some function c^, then we have 


02 — - j- — CgOl — 

TTjTT-' 


2Prl 
" jAWn 


Heia'^ 


(2.31) 


The second condition in (2.21), concerning the behavior of <52 will then be satisfied provided 
that 


d — £ 1+82 — s>p> 0 . 


(2.32) 


This differs from (2.25) by the last term in the left hand side of the inequality, which is 
defined as 


s = 


din cl 
dN 


(2.33) 


Although ( ^.30 ) is not completely general, it does cover some important cases. For instance, 
if the Lagrangian is linear in then the metric Pjj = Gij{(f>^) is just a function of the 
fields, and ( 2.30|) is satisfied with the speed of sound = 1. In this case di = S 2 and the 
and (2.32|) coincide. Another interesting example where (2.30) is satisfied 


conditions 

with a non-trivial speed of sound 7 ^ 1 is the case of a multi-field DBI action, which is 
Finally, Eq. 


discussed in |^] Finally, Eq. (|2.3C| ) also applies to generic one field models. It is easy to 
check that in such case, the attractor condition ( p.32 ) is related to the absence of a growing 
mode for the curvature perturbation on uniform field hypersurfaces 

Eor a generic multi-field model the condition that 62 vanishes in the asymptotic future 


will be more elaborate than ( 2.32 ), and should be worked out on a case by case basis. We 
leave this as a subject for further research. 


®More generally, when Pkl.ij is given by a linear combination of Pjj as 
Pkl,ij = JiPijPkl - I - f2{PiKPjL + PilPjk) , 

where /i and A are functionals of 4>^ and , AP and Ap are linearly related as in Eq. ( [2.3Cl| ). In such a 
case, the second attractor condition can be formulated as (2.32). 
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3. Separate universe approximation and 6 N formalism 

The approximate homogeneity and isotropy of cosmological evolution entails a useful re¬ 
lation between the curvature perturbation and the differential e-folding number, which is 
used in the so-called 5N formalism. Let us start by reviewing the conditions for the reali¬ 
sation of a universe with such approximate symmetries. The superpotential formalism will 
be useful in clarifying the conditions for the validity of this approach, particularly in the 
implementation of the momentum constraint. 


3.1 Separate universe approximation 

The 6 N method is based on the “separate universe assumption”. This is the statement 
that when a characteristic (physical) scale L of fluctuations is much bigger than the Hub¬ 
ble length H~^, i.e. L S> each region of the universe of Hubble size evolves as a 

locally homogeneous and isotropic FRW universe. This is justified if each Hubble patch is 
approximately homogeneous and isotropic up to corrections of order e <C 1, with 

1 

Then, since different Hubble patches are causally disconnected for local theories, they 
should evolve independently of one another. 

Consider the ADM line element, 

ds^ = —a^dt^ -|- 7jj(dx* -t- /3Mt)(dx'^ -|- , (3.1) 

where a and /3* are the lapse function and the shift vector, and 'jij is the spatial metric. 
Here, the shift vector /3* with the index i should be distinguished from the beta function 
with the index I. We parametrize 'yij as 

j.j = , tr = 0 . (3.2) 


The time-like congruence orthogonal to t = const, slices has a unit tangent vector given 
by —/?*), and its expansion K is given by 


K = = 


^9 


= a-^ d{H + 7^) - A/3 


(3.3) 


where Di is the covariant derivative defined with respect to the spatial metric The 
e-folding number may be defined as the integral of the expansion along the normal congru¬ 
ence; ^ 


Af — - J K a dt. 

^Alternatively, we may define the co-moving e-folding number AL = J{Kc/d)dT. Here dr = 
\J— PiP^dt is the element of proper time along a co-moving worldline, and Kc is the expansion of the co¬ 
moving congruence, whose unit tangent vector is given by n/f = 0)- T is straightforward 

to check that for dip^ = 0(e^), Me ~ M, up to terms quadratic in the spatial gradients. 



- 9 - 



In what follows, the metric is assumed to take the FRW form in the long wavelength 
limit e —>■ 0, and we shall use the gauge where the traceless matrix h is also transverse, 


d^hij = 0 . 


(3.5) 


For scalar fields, the anisotropic stress Tij — {1 /d)'^^^Tki'yij is of order O(e^). As shown in 
Ref. [551, using Einstein’s equations, after neglecting terms of lower order in e, which decay 
like a~'^, we have 

9# = 0(e2), (3.6) 


and 


kj = 0 {e^). 


(3.7) 


Using these conditions and introducing the derivative with respect to the number of e- 
foldings A7, 


Om ^ ^a, , 


the Einstein equations and the field equations of the scalar fields read |5(;] 

2 k^ 


K'^ = 


d{d-l) 


dMK = - 


K 


d-l 


- d^P) + 0 {e^) , 


KPijdj^(j)^dM4>‘^ + 0{e^) , 


Kdu {PijKd^(t>^) + dK^PijdM^^ - d^{dP/dcj)^) = 0{e^) 


(3.8) 


(3.9) 

(3.10) 

(3.11) 


At the leading order in the gradient expansion, these equations coincide with the back¬ 
ground field equations, where K should be understood as dPl for this comparison. Eor 
definiteness, in what follows we shall refer to Eqs. (|3.9|) - (|3.11 ), together with Eqs. ( |3.6| ) 
and (^), as the separate universe approximation. 

3.2 Momentum constraint 

The momentum constraint has no counterpart in the background field equations, and so it 
might enforce additional requirements for the validity of the separate universe approxima¬ 
tion. The momentum constraint is given by [p3] 


diK = — 


d-l 


KPudMck^dicj)^ + 0{ae^) . 


(3.12) 


A factor of a is inserted in O(ae^), since here we are considering the comoving derivative 
of the extrinsic curvature, while in our conventions a factor of e corresponds to a physical 
gradient. On the other hand, taking the spatial derivative of the Hamiltonian constraint 
and using Eqs. (^.lOD and (|3.11j ), we obtain 


diK = - 


d-l 


KPijdN(j>^dict>-^ + Bi + Oiae^), 


(3.13) 
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where Bi is defined as 


Bi = 


k^K 


^ {d-l)dMln{e'^K) {PijdM(t>^) " Qm [PijdM4>^) • (3.14) 


Here we used Eq. ( |3.10| ) to rewrite the denominator on the right hand side. Comparing 
(3.12) and ( |3.13| ) we see that the consistency of the Hamiltonian and momentum constraint 
requires that 

a-^Bi = 0{e^). (3.15) 


Sugiyama, Komatsu and Futamase pointed out that the condition ( 3.15 ) is automati¬ 
cally satisfied under the slow-roll approximation. Here, we argue that this conclusion is not 
restricted to slow-roll, but follows more generally from the attractor behaviour discussed 
in Subsection Indeed, repeating the same argument as in the background, we can 


express the field equations at the leading order of the gradient expansion with the use of 
the superpotential as 


^ ^ d-ldlnW 

PljdMft) — - — -jtIT-^ ^(^ ) ’ 




dcj)^ 


(3.16) 


where the superpotential is now related to the extrinsic curvature by K = 2dW + 0(e'^). In 
the attractor regime, the dependence of the superpotential on the integration 

constants ck can be neglected. Then, substituting Eq. (|3.16| ) in Eq. ( ]3.14|) , the leading 


terms in the two expressions within round brackets in the left hand side of Eq. (3.14) cancel 
each other, and we are left with a~^Bi = 0{e^). 

Beyond the attractor regime, we need to consider the dependence of W in the constants 
Ck- In that case, substituting Eq. (3.16) in Eq. (3.14) we have 


a ^Bi = - 


2dW 


dj\fln{e'^PW)dJ\f\ dcj 


a In IT 


(a '^diCj) + 


(3.17) 


where the total derivative with respect to M is taken along the dynamical trajectories at 
fixed Ck- From Eq. ( |2.1£ ), 


dW 

dcj 


= A 


J^-dAf 


+ Oie^), 


and after some simple algebra, we have 

a~^Bi = 2da~^A 


— Id _ „—d aJ , /n/',3\ 


+ 0(e"). 


(3.18) 


(3.19) 


The first term in the right hand side contains only one spatial derivative, and so it is naively 
of order e in the gradient expansion, rather than e^. Thus, unless the ck are constant in 
space, it might seem that the momentum constraint (|3.12) is inconsistent with the spatial 


derivative of the Hamiltonian constraint given in Eq. (3.12)®. Note, however, that the 


®For instance, it was stated in 0 that in single-field models the integration constant in the solution 
of the H-J equation should be constant in space for consistency between the Hamiltonian and momentum 
constraint. But as we argue here, this is not necessarily the case. 
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term of order e is accompanied by a decaying function which scales as® a while the 
terms of order e® include gradients of non-decaying contributions. In particular, for spatial 


dimension d > 2, the first term in (3.19D falls off with physical wavelength faster than 


and hence it can be safely ignored for present purposes. For any given co-moving scale, 
the initial conditions for the long wavelength evolution are generated at horizon crossing, 
and the first term in ( |3.19| ) will be negligible soon after that. 

As noted in Ref. [^] , the momentum constraint can always be satisfied in the gradient 
expansion by modifying (p.6|), so that /3* includes terms of lower order in e, starting at 
a/3® = 0(e“^). As mentioned before Eq. ( |3.6[ ), such lower order terms can be shown to 
decay as a~'^ for matter whose anisotropic stress is of order e^. Using the Hamiltonian 
and momentum constraints at the linearized order in the flat slicing we can obtain the 
relation |M] ^ 


dA/3^ + B, = 0, 


(3.20) 


where we used that 


dt= 


(3.21) 


which follows from Eqs. 


and (|3.8|) , and assuming a flat gauge where TZ = 0. Using 


Eq. O), we can relate the term of order e ^ in the expansion of the shift vector /3® to 


the hrst term on the right hand side of (|3.19|): 


9^/3® =-^AIUo + 0(e' 


(3.22) 


where AIUo = A'^Acj is a slowly varying function of position. After this mode decays, 
a~^Bi is formally of order e® and difi^ = O(e^), consistent with the separate universe 
approximation 


3.3 5N formalism 

In the separate universe approximation, the last term in square brackets of Eq. 
ishes as 

A/3‘ = O (e^) 


van- 


(3.23) 


in the long wavelength limit. Inserting Eq. 
we obtain 


into Eq. (^) and neglecting 0{e^) terms 


SN{t 2 , ti;x) = Af{t 2 , ti; x) - N{t 2 , h) ~ 7 l{t2, x) - x ), (3.24) 


where N is the e-folding number for the unperturbed background evolution. 


°For the case of canonical fields, such scale dependence of the leading term in the gradient expansion of 
Bi was first derived in Ref. js^ . 

^When we include the second term in Eq. (3.21), our equation (3.20) differs from the equation (49) in 

Ref. S. 
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With the help of Eq. ( 3.24 ), one can map the spatial distribution of the scalar fields near 
the time when the relevant scale crosses the horizon, S(pl{x), to the curvature perturbation 
at some final space-like hypersurface Eg near the end of inflation, TZ{tf,,x). The spatial 
distribution of the scalar fields should be specified on the “initial” space-like hypersurface 
S* where the spatial curvature vanishes. 


TZ{U,x) = 0. 


(3.25) 


For the final hypersurface Eg we have some freedom. For instance, we may take it to be a 
hyper surface of constant energy density. 


p{tp^,x) = pe = const. 


(3.26) 


Note that Eqs. (p.25) and (3.26) dehne the times t^{x) and tp^{x) as implicit functions of 


position. Alternatively, by a suitable choice of coordinates the hypersurfaces E* and Eg 
can be made to coincide with initial and final t = const, slices. The adiabatic curvature 
perturbation is defined as: 


at 


pe? 


X 


= 7 ^(^ 


Pe ) 


X) 


(3.27) 


Using ti = t* and t 2 = tp^ in Eq. ( p.24 ), from ( |3.25| ) and (3.27) we obtain the familiar 
relation between and the differential e-folding number 


atpe^x) = 6N{tp^,U;x) = 6N{tp^-,5(j)i{x)). 


(3.28) 


For the last equality we used the separate universe approximation and assumed that the 
inhomogeneity in the e-folding number can be determined from the initial distribution of 
the fields = (5(/>'^(t*, x). More precisely, the e-folding number from the initial to 

the final hypersurfaces should be computed by solving the Friedmann equation in each 
Hubble patch. In principle, this requires as initial conditions, so the initial time 

derivative of the field distribution is also needed. However, as discussed in the previous 
section, when the trajectory is an attractor, the dependence on the initial time derivative 
dies off rapidly and then the initial distribution can be expressed only in terms of a*. For 
the one held case, this is discussed in detail in the following subsection. 

If the inhationary trajectories in held space converge to a unique one i0 , then the 
adiabatic curvature perturbation C becomes subsequently constant in time, i.e. indepen¬ 
dent of the value of the density pe which dehnes the hnal hypersurface. However, in general, 
we also need to deal with entropy modes at the hnal hypersurface. These can be calculated 
along similar lines as C- 

In particular, choosing the uniform held slicing with = (j)^ = const, as the hnal 
hypersurface, which we may denote by one obtains a set of D independent gauge 


invariant variables which represent the curvature perturbation in the different 
slicings (see for instance |]T, ]^|): 


= 0 


, x) = n{j:^i,x) = ^p,x) = 5N^^\(i>^-5(i>i{x)). (3.29) 
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(3.30) 


Expanding x) in terms of the scalar field fluctuations 5(f){{x), we obtain 

OO ^ 

X) = Y, 

with 


n=l 


N 


in - 

J'l ...Jn 




T T • (3-31) 

Such partial derivatives of can be determined from the solution of the equation of 
motion for the local FRW universe as a function of initial field values in the vicinity of 
the background solution. 

The relative entropy perturbation is dehned as the difference between two 
(see e.g. Ill, 


x) = d x) — x) 




{J)( 


(3.32) 


which can of course be expanded in powers of 5cj)l^ by substituting Eq. (|3.30D . 

Finally, we note that the uniform Hubble (or energy density) slicing is given by 


D 


6 H ^dlnW^r K 
0 = ^ = 2 ^ =- 


D 


I=l 


dcj)^ 


d-1 




(3.33) 


1=1 


On the other hand, at linear order, the curvature perturbation TZ changes under the time 
shift t t + 5t SiS TZ —>• TZ — HSt, while the field fluctuations S(j)^ change as S4>^ —)• 
Scf)^ — {(5^/K)HSt. Using these transformations and Eqs. (|3.27) and (3.29), we have 


D 




D 






S(p^ 

'■W 


(3.34) 


1=1 ' 1=1 

Hence, at linear order we can relate the curvature perturbations C, and by means of a 
simple expression involving the (5 functions. 


3.4 Linearized perturbations in one field models 

In Eq. ( p.28|) we have neglected the dependence of 5N on the initial momenta tt^. Let us 
show more explicitly how this approximation is justified in the simple example of a one 
field model. First, let us show that the full linearized solution is recovered from the SN 
computation. The e-folding number can be written as 

r<t’e 

N= ^d(^. (3.35) 

J<t>* 4> 

In this subsection, where we focus on one field models, we omit the index for cj)^ and 
VT/. The variation SN is due to the variation of the initial value of the field and to 
the variation of the initial momentum dvr*. The latter corresponds to a variation of the 
integration constant Ac^ in the solution of the H-J equation. Thus, we have 

SN = + SN 2 , (3.36) 

(pufi 
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where the first term is constant, and corresponds to the well known constant solution for 
C on super horizon scales, while the second term is given by 


6N2 = Ac^ 


f 

J 6* 


d /H 


Hdc, 


dN. 


This can be cast in the form 


dN2 = 


AIT 

~W~ 


- 52 dN, 


where 62 is defined in ( ^.311) . Using O) , we have 

ATT 2 


W 


Ha<^ 


AlUo. 


(3.37) 


(3.38) 


(3.39) 


Now, by keeping track of the last term in ( |3.3D , which is decaying and has been neglected 
in the relation (|3.28) between and 5N, we get 


C = 5N + 


I 


H 


dN = 5N - 


ATT 


W 


-dN. 


(3.40) 


In the last step we have used Eq. (3.22), and we have kept only the term of order in the 
gradient expansion, neglecting the terms of order e^. Combining with Eq. ( 3.37 ) we obtain 


/ 


(2 = - 52dN = -2dAIUc 


/ 


Hsia^ 


dN. 


(3.41) 


This coincides with the long wavelength solution of the standard linearized equation of 
motion for perturbations ||5^, as expected. If the background is an attractor, then <52 is 
exponentially decaying with N, and so by choosing our initial time appropriately, C2 can be 
neglected altogether. This is the standard decaying mode of the curvature perturbation. 
In the single-field case, this confirms that we can neglect the dependence of 5N on the 
initial momentum vr*. 


4. Primordial spectra from superpotential 


In this section we consider the spectra of the curvature and entropy perturbations, respec¬ 
tively C and at t = te in the case of a separable product superpotential. Let us start 
by considering the more elementary cross spectra for the which are dehned in Fourier 
space as 

= {2T:Y5{k + k')V^^^\h). (4.1) 

Note that is related to the power spectrum of 5(f)^ at the time t*, 

{5<l^i{k)5ct>i{k')) ^ {2T:Y5{k + k')Vli{k), (4.2) 


by 


vf\k) = Vf^^(fc) , 


(4.3) 
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where we have used Eq. ( 3.30|) . 

Likewise, we can consider higher order correlation functions. The bispectrum of ^ is 
defined by 

= {27ry6{ki + k2 + k3)Bl^^^\ki, k2, ks). (4.4) 

This can be decomposed into two contributions, 

k2, k,) = Bl^i^\ku k2, ks) + B^^ifjku k2. ks) . (4.5) 

The first one corresponds to the intrinsic non-Gaussianity generated until around t = t* 
and is related to the 3-point function of the field perturbations at horizon crossing by 


{27Tfd{k^ + k2 + ks)B^l^^^\ku k2, ks) = N^j[^NyjN^jf\dcPi^{k,)6cpi^ik2)dcPi^{ks )). 

(4.6) 

The second one is generated by the super-horizon nonlinear evolution. Using Wick’s theo¬ 
rem, it is given by 

^{IJK) l,^\ _ ]\AJ) 


B 


^2, k,) = N^j’NyjN);]Vil^:{k,)Vil^J{k 2 ) + (2 perms), (4.7) 


where we have used the symmetry of with respect to the lower indices 


4.1 Separable product superpotential 


Let us focus on the case where the Lagrangian in Eq. (^) is such that the metric Pjj is 
diagonal and each element I only depends on the kinetic term and on the field 

P,J = 5uKiiX^,cj)^) , . (4.8) 

Moreover, we will assume that the superpotential is given by a separable product, i.e., 

D 

Wi(p^) = (4.9) 

1=1 


By the assumptions (^) and (^), Eq. ( p.l2 ) becomes 


Kl{(j)P'Kl) 


dN 


dcj)^ 


(4.10) 


where in the argument of Kj we have traded the field velocities cj)^ by their expression in 
terms of fields and conjugate momenta, through Eq. ( p.lOf ) (which is also separable in this 
case). Note that the momentum tt/ oc dW/dcf)^ will actually be a function of all fields . 
Hence, for the purpose of making Eq. ([4.1(]| ) separable, we shall further restrict to the two 
following cases. 

The first case is when the kinetic term is linear in X^^ so that Ki = Ki{(f)^). In that 
case, we can always use a field redefinition to bring it to the canonical form 

D 

P{X^,ct>^) = Y,X^ (4.11) 

i=i 
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The second case of interest is the one field case, D = 1, in which case we can use Eq. ( 2.12| ) 
to express vr^ as a function of 0, for an arbitrary P{X,4>). 

In these two cases, solving for the velocity and using the definition of beta function 
o, one obtains 




(4.12) 


Hence the evolution of 4>^ as a function of e-folding number N can be determined without 
being affected by the other scalar fields. The e-folding number between the slicing 6 (j)^ = 0 
at t = tf. and the flat slicing at t = is then given by 




/ Hdt = 

'i* J(f) 


r- K( 
^ 




m 


(4.13) 


If the trajectory is an attractor, then after the decaying mode becomes negligibly small, 
the variation in the e-folding number is given by 


= - 


Pi 


(4.14) 


with From this expression we can calculate the derivatives of with 

respect to the variations in Spi: 




= 


dfl 




* = 




<^Pi 


{PlfdN, 


<5 


(4.15) 

(4.16) 


and so on. Hence, the trajectory is constructed out of D independent trajectories, and 
there is no mixing among the different field components. 

It is known that 6 N can be analytically solved also in different examples. Under 
the assumption of slow-roll, it has been shown that the number of e-foldings N can be 
computed analytically when the potential is a separable product of potentials where 

each one depends on a single held (5^, or when it is a separable sum |5^. This has 
been used to compute the non-Gaussianity in multi-held models of inhation, in the case 
of slow-roll evolution in two-held |^] and in multi-held models [l60| , ^]. The approach of 
||5^ , |5^ has been extended beyond slow-roll, exploiting the H-J formalism in the case of a 
sum separable Hubble parameter |^, |^, Since by Eq. ( 2.11| ) the Hubble parameter 
is nothing but the superpotential W{cj)^) (up to a factor 2), this treatment can be 
straightforwardly applied to the case of a sum separable superpotential. The approach 
of [^] was applied also to the separable product superpotential (Hubble) case by Safhn 
in Ref. |^|. As discussed in Appendix P, where they overlap, our results agree with 
this reference. Note that as we discussed above, for the separable product superpotential, 
when the background is an attractor we can compute 5N without introducing integration 
constants along the trajectory, as in Ref. [^]. 
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4.2 Primordial spectra 

In what follows, we shall assume that the perturbations in the different fields 6 cj)l are 
uncorrelated at the time of horizon crossing: 


vyAk) = 6^-^v^jjk). 


I.J-1 


(4.17) 


In the multi-field case, Eq. (4.17) will hold provided that the linearized equations of motion 
for perturbations decouple from each other. For instance, when are exponentials, it 
is easy to show that we can go to a basis in field space where perturbations are decoupled 
from each other. The reason is that, as we shall see, we can always change variables so that 
the corresponding potential V{4>^) depends only on one of the fields, while the other ones 
are massless. In this example, the slow roll parameter ei need not be small. Note that for 
wavelengths well within the horizon, the field dependence in the potential is unimportant 
and perturbations of the different fields are effectively decoupled from each other. This 
suggests that in a more general setting the fields will be uncorrelated near the time of 
horizon crossing provided that £2^1- 

Using Eq. (4.17) in the power spectrum of Eq. (4.3), we have 




(4.18) 


Using Eq. ( 4.15 ), we obtain 

(fe) = (fe) , (A:) ^ {k ), 

For canonical scalar fields in d = 3 dimensions, the field spectrum V^j is given by 


(4.19) 


(4.20) 


and is the same for all I. 


The result ( 4.IS ) is also valid for the case of a single field with arbitrary speed of sound. 
In this single field case 


'P4>Ak) 




2k^Cs*Px* 


Px 

Px + 2XPxx ’ 


(4.21) 


where Cs is the speed of propagation of fluctuations and a * denotes sound-horizon crossing, 

k — 

Defining the cross spectra for the curvature and entropy perturbations by 


{ak)ak')) = i27ry5{k + k')V^{k), 

(4.22) 

{C{k)S^-^{k')) = {2TTf6{k + k')P^s^j{k ), 

(4.23) 

{S^-^{k)S^^{k')) = {2Trf6{k + k')VsiJsKL{k). 

(4.24) 


- 18 - 










and using Eqs. ( 3.32 ), ( 3.34 ) and (|4.19f ), we obtain 


n(t) = E 




1=1 




im 




Vf^s^jik) — d 


Ple/de K 


/djePe 


/3e^ {/3iy PI {PiY 


V,. 


VsiJsKL{k) = (f 


D 2 


i{dlM— djM){6KM — dLM)'P4), 


(4.25) 

(4.26) 

(4.27) 


where an index e denotes a quantity evaluated at t = t^. Thus, the only non-vanishing 
components of VgusKL^k) are 


VsiJsij{P) = -VsusJiik) = d^ 


{Pif ' {Pi) 


+ 


V. 


(4.28) 


Since the power spectrum of does not vary after the Hubble crossing, the auto-correlation 
of does not vary either. 

From Eq. ( 3.34 ), the power spectrum of C, is given by the sum of the conserved power 
spectra of C,^ as 

V(^{k) = ^RjV^^\k), . (4.29) 

Notice that when the trajectory converges at t = te, say in the direction of / = 1, satisfying 

(4.30) 


Pl:>Pi\ 


Pi 

for / 7 ^ 1, the amplitude of C, is determined solely by the one for / = 1 at t = t* as in the 
single field case as 


V^{k)-Vf\k)^ 


(4.31) 


Then, the influence of the components 7 7 ^ 1 does not explicitly show up in the power 
spectrum of C; while the change of these components can be still traced through their 
contributions to 77*. This property was pointed out by Garcia-Bellido and Wands in the 
different separable example, where the slow-roll approximation is employed p 8 [ . 

Similarly, the bispectrum generated by the super-horizon evolution in Eq. (^) can be 
rewritten using Eq. (|4.17 ), 


^ 2 , k^) = + (2 perms) 


Inserting Eqs. ([4.IS] ) and (|4.16|) into Eq. ([4.32|) , we obtain 


B 


tifjk„k„ks) = 6^^6 


_ xlJxlKdlnPi^ji), 


’ {k^V^'p {k 2 ) + (2 perms) 


dW 


(4.32) 


(4.33) 


The bispectrum for the entropy perturbations is trivially obtained from Eq. ( 4.33| ), since 
is linear in In principle, the bispectrum of C, can also be found from Eq. ( 4.33| ) 
by using the nonlinear relation between Q and 
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5. 6N and holographic inflation 


In Refs. [24, 25, the primordial spectra were computed holographically by means of the 
dual quantum field theory which lives on the three-dimensional boundary. The computation 
from holography may provide an alternative way to address the primordial perturbations 
generated during inflation. In this section, we reexamine the primordial spectra derived 
from the 6 N formalism, comparing them to the prediction from holography. 


5.1 Inflation from holography 

In this subsection, we briefly overview the way to compute the primordial perturbations 
from the dual boundary theory, following Refs. |^, |^, ^]. During inflation, spacetime 

is quasi-de Sitter, i.e., the de Sitter symmetry 50(1,4) is slightly broken by the time evolv¬ 
ing inflaton field. In order to provide the boundary QFT which is dual to the inflationary 
spacetime, we need to slightly break the conformal symmetry in R^, which is also 50(1,4). 
In particular, to address a model with D scalar fields, we consider a boundary QFT whose 
action is given by 


D 

5'qft[x] = 5'cft[x] + X] / ^^dg^Oi{x) , (5.1) 

1=1 


where dO^ is the d-dimensional invariant volume and y is the boundary field. The second 
term describes the deviation from the conformal field theory. Here, Oj{x) is a composite 
operator of y and are the coupling constants. Solving the renormalization group flow, 
we can compute the beta function. 



dlu/u ’ 


(5.2) 


as a functional of g^, where g is the renormalization scale. In general, since the different 
components I can couple with each other, solving Eq. (^]^) analytically is hard. However, 
in the case where the beta function is separable as (3g = [3g{g^), i.e., in the case where 
there is no correlation between the different components of , we can analytically solve 
Eq. (5.2) as in a QFT with a single deformation operator. In such a case, we can also 
compute the auto-correlation functions of 0 {x) as a function of g, analytically. 

Assuming that the wave functions of the curvature and entropy perturbations are 
related to the generating functional of the dual quantum fields as 


V’bulk 


J' 


= A Z, 


QFT 


c, 


J' 


(5.3) 


where A is a normalization constant, we can compute the correlators of C, and . Here, 
to include only independent degrees of freedom, we introduced the entropy perturbations 
s^' as s^' = 5^^^ with I' = 2, ■■■, D. In the boundary QFT, the primordial perturbations 
C and s^' should be treated as external fields. Once the wave function V’buik is specified, 
we can compute all the correlators for and s^'. For instance, the n-point function for 
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is given by 


r ^ 

{C{XI)C{X2) ■ ■ ■ C{Xn)) = DC Y\ DS^' IV'bulkl^ C(®l)C(®2) ' ' ' C(®n) • (5.4) 

I'=2 

The cosmological evolution in the bulk is described by the D scalar fields as a function 
of time and of the spatial coordinates. To describe the bulk evolution by means of the dual 
boundary QFT, these bulk quantities should be related to the quantities in the boundary 
QFT. One may expect that the time evolution in the bulk will be described by the RG 
flow in the boundary and then the time evolution of cj)^ will be determined by the RG flow 
of . When we specify the relation between the time in the bulk and the renormalization 
scale in the boundary, t = f(/r), and also the relation between cj)^ and , 


9 ^ = 9 ^ 


(5.5) 


with the use of Eq. 


the correlators of C and can be described by g^ and the 


correlators of O |25, 261. Then, the correlators of the primordial perturbations can be 
holographically computed by solving the RG flow in the boundary. 


In Refs. 


it was argued that in the de Sitter limit the renormalization 


scale g should be proportional to the scale factor a. 


g (X a. 


(5.6) 


However, the relations suggested in these references differ from each other when the solu¬ 
tions deviate from de Sitter spacetime. In Ref. [^], considering the RG flow with two fixed 
points (a fixed point (FP) is a point where the beta function vanishes), which corresponds 
to the time evolution in cosmology from one de Sitter to another de Sitter, it was shown 
that with the choice of Eq. ( |5.6| ), the power spectrum of the curvature perturbation C in 
single field models is conserved at large scales so that the holographic computation gives a 
result consistent with the standard cosmological perturbation theory. Meanwhile, a more 
subtle issue is left unresolved for the conservation of the bispectrum |^|. 


5.2 Comparison of the bulk and boundary computations 

In this subsection, we compare the result from the boundary computation, which is obtained 
by solving the RG flow, to the one from the bulk computation which is obtained in the 
6 N formalism. When the RG flow with D deformation operators is separable, i.e., it is 
given by the D copies of the RG flow with single deformation operator, using the conformal 
perturbation theory, we can derive the beta function f3g as [|^] 


I ^ dg^jg) 
^ dln/i 


„d/2 

(A^ - d)g\g) + ^ ’ 


(5.7) 


where is the scaling dimension of , c is the central charge, and Cj is the structure 
constant. Solving Eq. (^), we can compute g\g)- Once g^ = g\(t>'^) and t = t{g) are 
determined, g\g) gives the time evolution of the scalar fields (j)^ in the bulk. 
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Assuming the relation = ncj)^ and Eq. (|5.6|) leads to fig = (3 , where fi^ is defined 
in Eq. ( |2.17 ). In this case, one can compute the primordial power spectra of from the 
boundary QFT whose RG flow is separable. Under these assumptions, one finds 






WM 




(5.8) 


where V^i{k) is the power spectrum of in the flat gauge, computed from the two-point 
function of the boundary operator . Since fig = fi^, this result agrees with Eq. ( 4.ig|) . 

Moreover, in holography the power spectra of C and are related to those of (f^s by 
Eqs. ( 4.25| )-(4.27) with fig = fi^. Hence, the power spectra of C and computed from 
the boundary QFT agree with those computed from cosmological perturbation theory in 
the bulk. Notice that, in the case of separable trajectories, the conserved power spectrum 
of can be computed in the same way as in single field inflation. Therefore, the fact that 
both calculations agree on the conservation of {k) directly follows from the fact that 


they agree in the single field case |24, ^]. 

Finally, we note that instead of using the simple relation , one may identify 

the couplings g^ in the boundary with the scalar fields in the bulk by more non-trivial 
relations g^ = g^{(f>'^). In this case, the relation between the spectra of and those of the 
boundary operators may become more complicated, due to the Jacobian \\dg^/dcf>-^\\. Yet, 
changing the identification g^ = g^{(f)'^) can be simply understood as a field redefinition. 


6. Case studies 


In this section, we compute the primordial spectra for D canonical scalar fields, using the 
formula derived in Sec. For our purposes, we consider a separable superpotential as a 
product of exponential superpotentials, i.e., 


W{^^) = Wo exp 


D 


E/'U) 


( 6 . 1 ) 


where Wq is constant. Since the superpotential W{(f>^) directly gives the Hubble parameter 
as H = 2IU, the summation of over all I = 1, ■ ■ ■ , D should increase in time so that H 
decreases in time. 

It is instructive to consider the case where this superpotential describes D scalar fields 
with canonical Lagrangian, 


D 


P{X, cf^)=X- V{cf^) , ^ E 


i=i 


In this case Eqs. and (|2.12|) become 

2k^ 


= 




( 6 . 2 ) 


(6.3) 
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and 


II 


2 {d-l) dW{(j)^ 


dcj)i ’ 


H = 2W{(j)^). 


(6.4) 


Using these equations we can compute the potential U((/>^), which is given by 


V {<!>') = V„ 


7=1 






[d^^J 

exp 


D 


7=1 


lrj.l\ 


_ 2d{d-l)Wi 
ho = -:5- 


(6.5) 


Note that the potential is not, in general, separable. The beta function f3^{c/)^) is given by 

d(j)^ d-1 df^ 


In terms of the slow-roll parameter ej is given by 

d-i/d/^y 

^ \d(j)^ J 

6.1 Constant j3^: Power-law inflation 

First, we consider the case where is linear in 


( 6 . 6 ) 


(6.7) 


( 6 . 8 ) 


where pj is a dimensionless constant parameter. Then, the potential V{(j)^) becomes a 
separable product of exponential potentials for each (j)^ 




D 


7=1 


d-1 




(6.9) 


7=1 


In the single field case, this is known as power-law inflation Using Eq. (6.8) in 

Eqs. (|6)6[ ) and s one finds that /3^ and the slow-roll parameters sj become constant. 


/3^ = PI 


r? 

J _ Pi 
^1 — 


d-1 


( 6 . 10 ) 


In this case, the coupling constant of the dual boundary theory with (3^ = jSg blows up 
both at the IR and UV limits except for the trivial case with pi = d. 

Solving Eq. (|6.6|), we can compute the evolution of cf)^ and H as 


K,(j)^{N) = K(j)i +pi{N - iV*) 
H{N) oc , 


( 6 . 11 ) 

( 6 . 12 ) 


®The case where the potential is given by a separable sum of exponential potentials is known as assisted 
inflation, see for instance 
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where is an integration constant and we remind the reader that ei = 


expected, integrating Eq. (6.12) in time we obtain the power law evolution 

a{t) oc . 


(6.13) 


Inflation requires ei <C 1. 

Using Eqs. (4.25)-(4.27) and assuming that all flelds have the same power spectrum 
the power spectra of ( and are given by 


(k) = 0 , 




p, 




ik) = 


d — 1 


+ 




(6.14) 

(6.15) 

(6.16) 


where we used /3‘^ = {d— l)ei. The power spectrum of C is given by the same expression as 
the one for the single field case and the amplitude is frozen after t = For scale invariant 
field fluctuations in three dimensions V^^{k) = Hl/{2k‘^) and the spectral index is given 
by 


n, - 1 = -- 


4ei 


1 - 2ei 

In this case, the bispectrum vanishes. 

6.2 Linear P^{4>^) 

Next, we consider the superpotential where also includes the quadratic term, 

1 


(6.17) 




(6.18) 


where pi and qj are constant parameters. In this case, the beta function is given 

by the linear function as 


13^ =pi + 2qiK(l)^ , 


(6.19) 


and V{4>^) is given by 


U(</.^) = Uo 


D 


1 - 


d{d - 1) 


'^{pi + 2qiK(l)^f 


i=i 


D 


exp< 


d- 




1=1 


( 6 . 20 ) 


For qj 7 ^ 0, the potential V{cj)^) is not a separable product and it is not easy to 
analytically solve the Klein-Gordon equations, which are not separable. However, since 
the superpotential is a separable product, Eq. (|6.6|) can be easily solved, which gives 


K<P^iN) = K^i + (^K^i + - i) . (6.21) 
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As expected, in the limit g/ —>• 0 we recover Eq. ( |6.11 ). 
beta function (3^, 


Using Eq. (|6.2lD , we obtain the 


^\N) = Pi 629^^-^.) . (6.22) 

At late times, the beta function blows up for g/ > 0 and approaches 0 for qj < 0. In the 
perspective of holography, where the late time in cosmology corresponds to the UV limit 
in the boundary QFT, the boundary theory dual to the latter case has the FP in the UV 
limit /i —)• oo. During slow-roll inflation, /3^ should be kept much smaller than 1, requiring 

Pi £29/(^-^*) < 1. (6.23) 


Using Eq. (6.21) in Eq. (6.7) we hnd 


£[[N) = e[^ , £{^ = {pi + 2qiKc/)i)‘^ , 


(6.24) 


using which we have // = [{d — l)e{ — p]]/[i{d — 1 )(? 7 ]. Using this in Eqs. (6.1) and (^), 
we can give the following expression for the Hubble parameter 


H{N) = 2Wq exp 


E 

L I 


P] - (rf- l)£i(^) 

4(d - l)qi 


(6.25) 


For arbitrary values of qi with qi 7 ^ 0, the Hubble parameter H decreases in time, which 
is simply because ei > 0 . 

Since P^ varies in time for qi / 0, the power spectrum of ( varies also after t = t^. 
Given that the scalar fields take values pi sd. t = (more precisely, for each separable 
trajectory of p\ the final time is specihed by a value of each held, pi), the power 
spectrum for p is given by the summation of the conserved power spectra for as in 


Eq. ( 4.29 ) with the ratio Rj given by 


Ri = 


{pi + 2qiKpi) 




[Ej=i(PJ + 


2 • 


(6.26) 


Similarly, using Eqs. (4.26) and (4.27), we can compute the cross-correlation between p 
and and the auto-correlation of . In this case, since the beta function varies in 
time, given in Eq. ( |4.33|) , takes a non-vanishing value. 

6.3 Quadratic P^{p^) 

Next, we consider f\p^) which includes a cubic term as 


1 


= JZrj + ri{Kp^)^] 


(6.27) 
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where pi, qj, and rj are constant parameters. Now, the beta function and the 

potential V{c/)^) are given by 


P )= PI + 2qiK(l) + 3ri{K(j) ) , 

D 


(6.28) 


V{<P^) = 1^0 E [P^ + 


D 


X exp< — 


d- 




1=1 


For later use, we introduce 


Dj = qj - 3rjpi. 


(6.30) 


For Dj > 0, the beta function /3^ vanishes at two different values of cj)^, i.e., the dual 
boundary theory has two FPs. For Dj = 0, /3^ vanishes at one value of (j)\ i.e., the 
boundary theory has one FP. For Dj < 0, the beta function does not vanish at any 
values of (j)\ i.e., the boundary theory has no FPs. Using the beta function, given in 
Eq. (6.28), we obtain the power spectrum of C as in Eq. (4.29) with 

[pi + 2qiK(j)i + 3ri{K4)lf]‘^ 


Ri = 


{Ej=i [pj + + 3rj(K0/)2]^| 


(6.31) 


In the following, we study the background evolution of these three cases in turn. 


6.3.1 Dj < 0: RG flow with no fixed point 

Eirst, we consider the case where does not vanish. Solving Eq. (^), we obtain 


K/(iV) = + ^^tan0^(iV) 

3ri drj 


with 


e^{N) = y^l{N-N,)+tan-^ , Ei = qi + 3riKcjii 


(6.32) 


(6.33) 


Inserting this solution into Eq. ( |6.28 ), we can compute the time evolution of the beta 
function (3^ as 


/3'(iV) = -^[l + tan2 0^(Ar)] 


(6.34) 


The beta function starts to grow rapidly, when tan0^(iV) becomes 0(1). Now, the Hubble 
parameter, given by 


H{N) = Ho exp 


L /=! 


^27(d-l)rf 


D 

E 


tan0^(iV) (tan^ 6\N) —3) 


(6.35) 
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with 


Hq = 2VFo exp 


■ D 

E 


L/=i 


qi 

27{d - l)rj 


{9rjpj - 2qj) 


(6.36) 


decreases monotonically in time. 

In this case, after the slow-roll time evolution, inflation ends when 9i{N) ~ 7r/4 for at 
least one of the Is and, afterwards, the Hubble parameter starts to decrease more rapidly. 
Therefore, this case can provide a graceful exit to inflation. Meanwhile, in the boundary 
side, the dual QFT does not have any FPs and the RG flow is dominated by irrelevant 
deformations in UV. It may be interesting to study such boundary QFT. 


6.3.2 Dj = 0: RG flow with one fixed point 

Next, we consider the case where the beta function (3^ vanishes only at Kcj)^ 
this case, the time evolution of (N) and (3^ (N) are given by 




^ _1_ 

3r/ 3ri{N - N,- EJ^)' 


and 


f3\N) 


1 1 
^(N-N,- Ejy ' 


-qi/3ri. In 


(6.37) 


(6.38) 


At late times with {N — — Ej^) ^ 1, c/)^ approaches the constant value —qj/3rj, where 

the beta function f3^ vanishes. In this case, the boundary theory has one FP in the UV. 
When all components satisfy Dj = 0, the Hubble parameter becomes constant at late 
times and the universe becomes the de Sitter spacetime. This solution does not provide a 
realistic model of inflation because there is no graceful exit. 


6.3.3 Di > 0: RG flow with two fixed points 

Finally, we consider the case where the beta function f3^ vanishes at two different values: 


I —qi ± 


K(t)^ = 


In this case, solving Eq. (|6.2j ), we obtain 


and 


with 


ct>\N) = 


I3\N) = - 


3r/ 

(Pi + 

1 + ’ 

AD I 

3ri [l + e2®'W]2 


(6.39) 


(6.40) 


(6.41) 


e\N) = ^l{N , Ei = qi + 3riKcPi. (6.42) 
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Notice that both in the early time and late time limits, where 9 ±oo, the beta function 
I3\N) vanishes, approaching the de Sitter spacetime. In this solution, takes the constant 
values (jiL and </>;![_ in the limits N —oo and N oo, respectively. If all components of 
the scalar fields satisfy Dj > 0, the solution describes the transition from one de Sitter to 
another de Sitter. In this case, the dual boundary theory has 2 FPs both in the IR and 
UV and its RG flow is driven only by relevant deformations. Such boundary theory was 
studied by means of the conformal perturbation theory in Refs. 

If all components satisfy Dj > 0, the universe becomes the de Sitter spacetime at late 
times. Only if there exists at least one component I with Dj < 0, inflation can terminate 
as discussed in Sec. Notice that the components with Dj > 0, whose decrease in 

time, will satisfy |/3/e| <C l/Sjgl at sufficiently late times. Then, Rjs become negligibly small 
for I ^ I and hence they do not explicitly contribute to the spectra of i^, while they can 
still contribute implicitly through the Hubble parameter and the slow-roll parameters at 
the Hubble crossing time. 

7. Conclusion 

In this paper, we reviewed the superpotential formalism for multifield inflation, and ex¬ 
tended it to include the case of non-minimal kinetic terms. The superpotential is useful in 
characterizing the attractor behaviour of inflationary trajectories, as well as to assess the 
validity of the separate universe approximation. Furthermore, the logarithm of the super¬ 
potential plays an interesting role in the dual description of inflation, as the c-function for 
the RG flow in the boundary theory, whose gradient is related to the beta functions. 

Using the 6N formalism, we obtain simple expressions for the power spectra for adia¬ 
batic and entropy perturbations in the case when the superpotential is given as a separable 
product. In that case, the trajectory for each field is convergent even when the whole 
trajectory is not, and the power spectra can easily be found by solving the corresponding 
separable background trajectories. 

The bulk solution for the separable product superpotential corresponds to a boundary 
QFT with a separable RG flow, where the deformation operators are mutually uncorrelated. 
In such case, we showed that the power spectra of the adiabatic and entropy perturbations 
computed from the 5N formalism agree with the ones computed by solving the RG flow of 
the dual boundary theory. 

The separable case we addressed in this paper is described by D copies of the single field 
case. The power spectra of the adiabatic and entropy perturbations can be expressed in 
terms of such single field power spectra by linear relations. Because of that, the agreement 
of the bulk and boundary computations for the curvature perturbation in the single field 
model directly implies the agreement for the adiabatic and entropy perturbations. It would 
be very interesting to check the agreement in more non-trivial multi field models. 

Finally, with a view to phenomenological applications, we have considered some case 
studies of RG flows with a polynomial c-function, with terms up to quadratic order in the 
fields. These contain a range of possible behaviours from the infrared to the UV, which 
may hopefully illustrate the results which should be expected in more realistic scenarios. 
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A. Second order equations of motion from the superpotential 

Here, we show that the usual second order equations of motion (2.5) and ( p.Tj ) follow 
from the first order equation of motion (2.12) with (2.11), for any superpotential W which 
satisfies the H-J equation ( |2.13|) . 

First, we square (2.12|) to obtain 


PijPklX^^ = 


Here, we are using the notation 


Wi = 


dW 
d4>^' 


(A.l) 


(A.2) 


In what follows, we assume that Pjj is an invertible matrix and we shall denote its inverse 
as P^'^. We use these matrices to raise and lower the field indices. In particular 


XiK — PuPriX-^^ = — — y^WiWk, 


and we define 


X = PijX^J = P^^Xij. 


IJ 


The H-J equation (2.13) can then be written as 


^-^W^ = 2X-P, 


(A.3) 

(A.4) 

(A.5) 


where it is understood that any occurrence of is replaced by its expression in terms of 
cj)'^ and Wj{(j)'^) through Eqs. ( |2.1C ) and (2.12). 

Taking the total derivative of (A..5) with respect to (j)^, and using (2.12), we have 


-j dX dP „ dA^-^ 

-dHPu<t. = 2^ - ^ 


(A.6) 


where the 4>^ dependence of P{<f>^, X^'^) has been separated into its explicit dependence 
and its dependence through field velocities contained in X^'^. Using 

^ dA^^ dA ^KLdPKL dA , ^ dP^^ 

d(t>^ dcj,^ ^ d4>^ d^^ + d4>^ ’ 
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and 


we have 


d(/>^ 

dP 


= P 


KL 


dXKL , 


+ 


d(f)^ 


dcj)^ 

= dHPij^-’ + P 


Xkl, 


dcp^ 

.KL dXxL 


We can evaluate the last term as follows: 

pKl^^KL ^ j,KL ^d-lf d-^W ^ 

tP iT 


d4>^ 

-2{d-l)dWi 
K? dt 


(A.8) 

(A.9) 


_d 

dt 


. (A.IO) 


dcj)^ k‘^ d(j)^d(j)^ 

With this, Eq. (A.£) coincides with the second order equation of motion for (j)\ Eq. (p.5|), 
while the Eriedmann equation (2.7) is also satisfied because W solves (A.5). 

B. Dilution of dW/dcx with cosmic expansion 

In this appendix, we find the time dependence of the derivative of the complete solution 
of the H-J equation, W{4>^ ,ck), with respect to the integration constants. Erom the 
Friedmann equation o , we have 


Ad{d-l)^^dW dX^^ 

= PKL- 


(B.l) 


K- dcj dcj ’ 

where we note that p only depends on ck through the kinetic variables X^^. Now, from 


(A.3), we have 


PkiPljX^^ = 


where, again, we are using the notation (A.2). Taking derivative of 
CK, and then contracting with the ’’inverse metric” P^'^, we immediately find 

dxLM 


(B.2) 

with respect to 


{^Pij,lmX^'^ + Plm) 




K 


(B.3) 


dcK tP dc 

Noting that p = 2PjjX^-^ — P, Eq. (|B.3D can be rewritten as 

dX^^ (d-l)ddW 
dcK dt^’ ^ 

where we used O) to express the right hand side as a total time derivative. Here plm 
denotes the symmetrized derivative of p with respect to X^^. 


Substituting Eq. (BH) in Eq. (^) we obtain 

BW 1 d dW 


and using H = 2IT, we have 


which leads to 


W 

dt 


dcj 2d dt dcj ’ 
dW" 


In 


dcj 


= -dH, 


dW 

dcj 




(B.5) 

(B.6) 

(B.7) 


where are integration constants. Hence, the dependence of W on the integration con¬ 
stants becomes smaller in time along the dynamical trajectories, diluting with volume. 
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C. Alternative way to compute 6N 


In this appendix, we derive the primordial spectra of ^ in an alternative way, following 
Refs. 


59 


When several scalar fields contribute to the background evolution, the 
background equations of motion can be solved analytically only for some special cases, e.g., 
the case where the background evolution can be described by separable equations |^, M 


^]. In Ref. 1^], Garcia-Bellido and Wands computed the power spectrum of C for the two- 


field model whose potential is given by a separable product as ^ 2 ) = Vi{4>i)V 

under the slow-roll assumption. In Ref. it was shown that a similar analysis can be 
done also for a separable summation potential (^ 2 ) = Rl(0i) + Vi4‘2)- This analysis 

was extended to the case with arbitrary number of the scalar fields @] (see also Ref. [^] ). 
In these discussions, it’s crucial that the Klein-Gordon equations for the scalar fields can 
be recast into the first order equation by employing the slow-roll approximation. 

With the use of the superpotential W{(j)^), which is related to the Hubble parameter 
as H{cj)^) = 2W((/)^), the field equations for the scalar fields can be recast into the first 
order equations without employing the slow-roll approximation. In Refs. [32, 


, under the 

assumption that the Hubble parameter is given by the separable summation as = 

Byrnes and Tasinato computed 6N, following the method by Garcia-Bellido 
and Wands |58|. This analysis was extended to a non-canonical scalar field in Ref. |]6^ 
and to the case with the separable summation superpotential (or the Hubble parameter) 
in Ref. |^. In Appendix, we summarize the computation of 5N in the method by Garcia- 
Bellido and Wands 


for the separable product W The overlapped part with Ref. |^| 

agrees. 

When the superpotential W{4>^) is given by the separable product as in Eq. (^) and 
P/j becomes diagonal as in Eq. (fl.Sj), the beta function (3^ is given by a functional of as 


(C.l) 


To solve the background trajectory, we introduce the integrals of motion C, as 


Ci = 


f dcj)^ f 

J (</>') J 


d(l)‘ 






,i+i^ 


(G.2) 


with i = 1, ■ ■ ■ , D — 1. Using Eq. (CT), we can verify that Ci actually stays constant along 
the trajectory. With the aid of the constant parameters Ci, the change of the e-folding 
along the trajectory can be expressed only by one of the fields, say as 


N 


{te, C, {Ci}^-,^) = j'j Hdt = 


K.d(l^ 

WW) 


(C.3) 


On the second equality, we used dt = d(j)^/cj)^ and Eqs. ( 2.11 ) and (CT). Unlike 6N^, 
which is determined only by 6(l)i for the separable product case, the e-folding number 
( C.3 ) depends also on (j)^ with / / 1, since (/>g depends on CjS, which are determined by 


using (1)1 with I = 1, ■ ■ ■ , D. Now, taking the variance with respect to (1)1, we obtain 


diV = - 


Kd(j)l 

PI ' A 


Af 


1 V—V , ,j 


(C.4) 


1=1 
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In the following, we compute dcj)l/d4>l. The integrals of motion Ci can be expressed 
in terms of (j)i and hence we obtain 


E 


dct>i da 


_ 

d<t>i ^ da d4>i 


(C.5) 


Using Eq. (C^), we obtain 


d4,i ssa^ 


(C,6) 


Choosing the uniform Hubble slicing, we specify the final time te as the time when the 
Hubble parameter takes a particular value as 


D 




(C.7) 


i=i 


Taking the derivative of with respect to Ci and dividing it by Hg, we obtain 


/=! 


dCi 


(C.8) 


Next, introducing 


i-i 


Ci = Y,C, 


2 = 1 


WW)~J WW) 


dcf)^ 


(C.9) 


with / = 1, • • • , I?, we compute d4>l,/dCi. Taking the derivative of Cj with respect to Ci, 
we obtain 


dCj 1 d(j)l 1 9(/>g 


(C.IO) 


where we noted that values of cj)^ depend on Ci chosen for each trajectory. Equation ( |C.10 ) 
is recast into 


d(j), 






dCi dc^ 


1 d(l)l 


^ -Qii 


with 


(C.ll) 


P) = ^ ^ j 1 (^ < 7 - 1) 

- dCi |o(i>/-l) 


(C.12) 


Using Eqs. ( |C.8D and ( |C.11D , we obtain 


d(pi _ oi 

dCi 


Ey=.+1 /3/e/3i 
PI 


- Qii 


(C.13) 
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with 


N 


i=i 


Using Eqs. (C^) and (|C.13 ), we obtain 

A 

d(l>i A 


1 r 

e_ 

I 


PlePi ^ 

~w - 


Inserting Eq. ( p.l5|) into Eq. (|C.4D , we arrive at the compact expression: 


D 


^ Plef^i 


1=1 


Al Ai 


Using Eq. ( C.16 ), we can obtain 


Ni = 


Nij = 


dN 


A/eAe 


d(t>l 

d'^N 

d(j)id(j)i 


Al Ai 


= Sij 


A/eAe dA^ 
PI (Ai^diV*’ 


(C.14) 


(C.15) 


(C.16) 


(C.17) 

(C.18) 


and so on. As in the case with the separable summation W{<f>^), can be immedi¬ 

ately computed for a given Wj{cj)^) and take non-vanishing values only if /i = • • • = /„. 
The power spectrum of the curvature perturbation computed from Eq. ( C.17 ) agrees with 
Eq. (14351) . 
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